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Abstract 

Background: Green microalgae represent a renewable natural source of vitamin E. Its most bioactive form is the 
naturally occurring /?/?/?-a-tocopherol which is biosynthesized in photosynthetic organisms as a single stereoisomer. 
It is noteworthy that the natural and synthetic a-tocopherols are different biomolecular entities. This article focuses 
on /?/?/?-a-tocopherol production in Stichococcus bacillaris strain siva201 1 biomass in a bioreactor culture with methyl 
jasmonate (MeJa) elicitor. Additionally, a nonlinear mathematical model was used to quantitatively scale-up and 
predict the biomass production in a 20 L balloon bioreactor with dual variables such as time and volume. 

Results: Approximately 0.6 mg/g dry weight (DW) of /?/?/?-a-tocopherol was enhanced in 5. bacillaris strain siva201 1 
biomass with the MeJa 50 liL/L for 24 hrs elicitations when compared to the control. The R 2 value from the 
nonlinear model was enhanced up to 95% when compared to the linear model which significantly improved the 
accuracy for estimating 5. bacillaris strain siva2011 biomass production in a balloon bioreactor. 

Conclusions: 5. bacillaris strain siva2011 is a new green microalga which biosynthesizes significant amounts of 
/?/?/?-a-tocopherol. Systematically validated dual variable empirical data should provide key insights to multivariable 
or fourth order modeling for algal biomass scale-up. This bioprocess engineering should provide valuable 
information for industrial production of /?/?/?-a-tocopherol from green cells. 
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Background 

a- tocopherol is a lipid soluble small molecule and 
the biologically active form of natural vitamin E. RRR-a- 
tocopherol is exclusively biosynthesized by photosyn- 
thetic organisms or green cells including algae, plants, 
and cyanobacteria [1-3]. Plant-based products are a pri- 
mary source of ^ftR-a-tocopherol in the human diet. 
For example, hazelnut is one of the richest sources of 
vitamin E [4]. It is known that vitamin E plays an im- 
portant role in human nutrition as a natural antioxidant. 
It was recently proposed that vitamin E is active against 
oxidative stress-related diseases [5]. Reportedly, it sup- 
presses telomerase activity in ovarian cancer cells [6]. 
Vitamin E enhances IL-2 production, gene expression, 
and is an effective therapeutic adjuvant [7]. Vitamin E 
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deficiency affects both T and B immune cell functions 
[8], the a-tocopherol transfer protein (a-TTP) gene, and 
neurologic dysfunctions. In animal models, vitamin E 
mixtures inhibit colon, lung, mammary, and prostate 
carcinogenesis [9], as well as prevent diabetes [10]. Cos- 
metic industries also extensively use vitamin E in skin 
care products. In addition, .ftftR-a-tocopherol prolongs the 
shelf life of meat [11]. Thus, the use of .ftftR-a-tocopherol 
continues to increase in nutraceuticals [12]. 

Bioreactor technology is the key component for the 
industrial scale production of bioactive small molecules 
for pharmaceutical applications. For instance, bioreactor 
technology was successfully developed and scaled-up to 
10,000 L for commercial-scale production of ginseng roots 
used for human health-related applications [13,14]. This 
reactor configuration has also been tested for RRR-a- 
tocopherol production in lab-scale photosynthetic hazel- 
nut root culture [15]. However, knowledge regarding the 
bioprocessing of green algal cells for production of RRR- 
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a-tocopherol in balloon bioreactors is limited. Previously, 
Stichococcus bacillaris strain siva2011 biomass was scaled- 
up in a lab-scale balloon bioreactor (4 L - 8 L), and a lin- 
ear fitting model for predicting scale-up was proposed 
[16]. The S. bacillaris strain siva2011 has unique lipid 
(Figure 1) [16] and vitamin E metabolisms to lead to bio- 
active .ftftR-a-tocopherol production. A nonlinear model 
enables more accurate estimates and should provide 
insight for quantitatively predicating algal biomass accu- 
mulation in a large volume bioreactor. The objectives of 
this study were to: 1) evaluate RRR-a- tocopherol content 
from S. bacillaris strain siva2011 biomass with Meja elicit- 
ation; 2) quantitatively predict S. bacillaris strain siva2011 
DW in a 20 L bioreactor with simplified stepwise dual var- 
iables (time and volume) from 4 L and 8 L data using non- 
linear regression. This systematic study could provide 
insight regarding stepwise nonlinear scale-up of S. bacil- 
laris strain siva2011 biomass for .ftftR-a-tocopherol pro- 
duction in a balloon bioreactor. 



Results and discussions 

Natural vitamin E occurs in two general forms, namely 
the tocopherols and tocotrienols, which are collectively 
called tocochromanols. Each has four distinct isoforms 
having a, |3, y, or 6 substitution on the chromanol. The 
natural a-tocopherol contains three methyl groups on 
the chromanol moiety at positions 5, 7, and 8. The phy- 
tyl or saturated side chain is attached to the C-2 position 
of the chromanol ring which has three chiral centers 
with a single RRR stereoisomer (Figure 2). The chroma- 
nol groups have two fused rings, a phenol and a tetrahy- 
dropyran, sharing a 2 carbon bridge [17]. These rings 
are moderately polar, giving them an affinity for the cel- 
lular membrane surface while the phytyl tail is hydro- 
phobic and normally associated with membrane lipids 
[18]. These structural features of &&R-a-tocopherol are 
efficiently acted on by the human hepatic a-TTP which 
is responsible for maintaining plasma a-tocopherol con- 
centrations [19]. 
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Figure 1 Gas chromatography mass spectrometry profile of fatty acid methyl esters from S. bacillaris strain siva2011, biomass from 
20 L bioreactor, working volume 8 L with 0.2% C0 2 on day 6. Methyl palmitate (CI 6:0), methyl hexadecatrienoate (C1 6:3), methyl oleate 
(CI 8:1), methyl linoleate (C18:2), and methyl linolenate (C18:3). 
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Figure 2 Chemical structure of /?/?/?-a-tocopherol. 
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The synthetic a-tocopherol called all-racemic-a-tocopherol 
is not identical to &£R-a-tocopherol. It is an equimolar 
mixture of eight stereoisomers which possess three 
chiral centers at positions % 4', and 8', giving rise to four 
diastereoisomeric pairs of enantiomers such as RRR, 
RSR, RRS, RSS, SRR, SSR, SRS, and SSS [20]. Moreover, 
a-TTP has a high affinity to a- tocopherol and has 
a 3-fold greater binding half-life when compared to syn- 
thetic a-tocopherol [21]. Thus, the bioactivity and the 
relative safety are different. Human proteins such as en- 
zymes and receptors usually exhibit high stereospecificity 
[20]. Therefore, the natural .ftftR-a-tocopherol is more bio- 
active than the synthetic form. 

i^^-a-tocopherol plays a major role as an antioxidant 
which prevents lipid peroxidation. In the photosynthetic 
cells, it may protect photosystem II during photoinhibi- 
tion and repair chloroplast mechanisms [18]. This is due 
to the hydroxyl group on the C-6 position which is the 
active site that donates a hydrogen atom. The phenolic 
hydrogen atom is capable of scavenging lipid peroxy radi- 
cals and quenching singlet oxygen [22]. RRR-a- tocopherol 
is recycled in the photosynthetic cell by cytosolic ascor- 
bate which oxidizes one-electron from the tocopheroxyl 
radical thus regenerating vitamin E [23]. This mechanism 
might protect the cell membranes. 

i^^-a-tocopherol is biosynthesized in photosynthetic 
cells via two different pathways [24], The phytyl domain 
precursor comes from an isoprenoid pathway, and the 
chromanol domain precursor comes from an alternative 
shikimate pathway homogentisic acid via complex en- 
zymatic reactions [25]. ^ftR-a-tocopherol is found in the 
chloroplast envelope, thylakoids, and plastoglobuli of the 
plastid. The vitamin E biosynthetic pathway has been 
elucidated in Arabidopsis [26]. The genes associated 
with vitamin E biosynthesis in photosynthetic organisms 
have been well described in literature [3,27]. A signifi- 
cant metabolic engineering effort has been made to im- 
prove vitamin E content both in plants [28] and in 



cyanobacteria [29,30]. Moreover, tocopherol production 
in plant green callus [31], cell [32], and root cultures 
[15] have also been reported. 

/?/?/?-a-tocopherol production in S. bacillaris strain 
siva201 1 

Photosynthetic algae are a potential alternative for pro- 
duction because they biosynthesize an abundance of 
i^^-a-tocopherol. For instance, freshwater Euglena gra- 
cilis [33], marine Dunaliella tertiolecta and Tetraselmis 
suecica [34], model system Chlamydomonas reinhardtii 
[35], and commercial algae Spirulina platensis [36] have 
been used to biosynthesize ^ftR-a-tocopherol. The new 
green alga, S. bacillaris strain siva2011 [16], produces 
significant amounts of ^ftR-a-tocopherol. This microalga 
has efficient photosynthetic mechanisms which facilitate 
the quick biosynthesis of vitamin E. 

Tocopherol production can be enhanced by molecular 
elicitation which is non-transgenic. Meja is a plant stress 
volatile signaling molecule which up-regulates several 
defense-related genes [37]. In plants, jasmonates are bio- 
synthesized via the octadecanoid pathway; exogenous 
Meja treatment is up-regulating secondary metabolic 
pathway genes especially those encoding for stress pro- 
tection [38]. Therefore, it is used as one of the potential 
molecular elicitors in plant root culture to enhance sev- 
eral pharmaceutical molecule productions [15]. For in- 
stance, Meja elicitation increases the activity of tyrosine 
aminotransferase in plant green cell culture which is one 
of the initial step enzymes involved in tocopherol bio- 
synthesis [31,39]. To increase &£R-a-tocopherol content 
in S. bacillaris strain siva2011, Meja elicitation was also 
used. The S. bacillaris strain siva2011 characteristics, cul- 
ture conditions, and bioreactor experimental designs were 
previously reported for lipid production [16]. 

Figure 3 illustrates the ^ftR-a-tocopherol content of S. 
bacillaris strain siva2011 biomass under elicitation with 
various concentrations of Meja. The unelicited culture 
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Figure 3 ffffff-a-tocopherol content in S. bacillaris strain siva2011, unelicited and methyl jasmonate elicited culture on day 4, 5, and 6. 



accumulated 0.7 mg/g DW of .ftftR-a-tocopherol which 
was detected starting during the early exponential growth 
phase. The lower concentration of MeJa, 50 (iL/L, after 
24 hrs elicitation, enhanced the production of RRR-a- 
tocopherol to the highest concentration, 1.3 mg/g DW. 
The higher MeJa concentrations or longer elicitation pe- 
riods were inhibiting both to the biomass growth and the 
resultant &&R-a-tocopherol production. MeJa can diffuse 
to cells either by intercellular migration or while in the 
vapor phase [40]. In plants, MeJa can transport from 
leaves to roots [41]. The vapor signaling can be trans- 
ported to distal plants via air, and the intercellular sig- 
naling can be transported via vascular process [42]. For 
instance, in an in vitro root culture the MeJa elicitation 
could trigger the defensive molecules accumulation via 
intercellular transport [13-15], whereas in in vivo plants 
during herbivore attack it can act as a volatile signal 
[42]. In the algal culture, MeJa elicitation could trigger 
the .ftftR-a-tocopherol production by either intercellular 
signaling or both. The optimum concentration can up- 
regulate the tocopherol biosynthetic pathway enzymes in 
S. bacillaris strain siva2011 which could increase the anti- 
oxidant. MeJa elicitation rapidly activates the defensive 
genes which also down regulates the photosynthetic sys- 
tem genes [38]. In addition, MeJa induces reactive oxygen 
species (ROS) which alter the mitochondrial and chloro- 
plast dynamics [43]. Thus, the higher MeJa concentrations 
or longer elicitation periods can produce uncontrolled 
ROS which can precede chloroplast or photosynthetic 
dysfunctions which could be inhibiting the tocopherol 
biosynthesis and biomass accumulation in S. bacillaris 
strain siva2011. In green cells, chloroplasts are an essential 
organelle for energy capture and transduction; a decline in 
photosynthetic activity is closely related to the decrease in 
the biomass. The typical MeJa elicited cells' symptoms 



were loss of chlorophyll, which causes the decline in 
the net photosynthetic rate, and degradation of 
ribulose bisphosphate carboxylase, etc. [44]. For in- 
stance, 100 \iL MeJa at 9 hrs elicitation altered chloroplast 
morphology and function which is associated with cell 
death [43] . Even though .ftftR-a-tocopherol indirectly reg- 
ulates the amounts of jasmonic acid [18], the decline in 
chloroplast efficiency could down regulate the RRR-a-toc- 
opherol metabolism. This suggests that higher MeJa con- 
centrations or longer elicitation could be cytotoxic beyond 
what was studied. 

Nonlinear regression 

Compared to green adventitious roots, photosynthetic algae 
have a higher ability to biosynthesize &&R-a-tocopherol in 
bioreactor cultures. The balloon bioreactor is a liquid-phase 
reactor with enhanced geometry and efficient fluid flow 
dynamics which could help provide higher mass transfer 
efficiency [14]. When compared to other photobioreac- 
tors, the balloon bioreactor had a larger headspace which 
efficiently captured light and enhanced photosynthesis 
[16]. Thus, S. bacillaris strain siva2011 was scaled- up in a 
balloon bioreactor (Figure 4) to investigate enhanced bio- 
mass accumulation. Of the four concentrations of C0 2 
tested, 0.2% yielded the highest biomass of 3.45 g/L in 4 L 
and 3.79 g/L (DW) in 8 L on day 6 [16]. The RRR-a-toc- 
opherol production was unchanged by the 0.2% C0 2 (data 
not shown). 

To quantitatively predict larger-scale algal biomass 
production based on lab-scale test data, the following 
stepwise structured approach was proposed [45,46]: 1) 
to set up the simplest model to linearize with dual vari- 
ables; 2) to model a nonlinear regression with dual vari- 
ables; and 3) to demonstrate multiple variables based 
on a nonlinear regression. Although scale-up predication 
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requires multiple variables, the 6 days algal culture does 
not significantly utilize all the media components. There- 
fore, selection of important dual variables can give insight 
on the efficiency of parameter selection for lab-scale valid- 
ation and initial scale-up prediction. In addition, the mul- 
tiple linear regression models might not incorporate the 
underlying nonlinear relationships [47]. So, in this study 
the second systematic approach was conducted to evalu- 
ate nonlinear modeling for scale-up prediction in a 20 L 



balloon reactor using dual variables. Nonlinear regres- 
sion models are more important tools than linear 
models because they provide better parsimony, inter- 
pretability, and prediction [48]. Figure 5 illustrates non- 
linear modeling of predicted S. bacillaris strain siva2011 
biomass accumulation in a 20 L balloon bioreactor with 
0.2% C0 2 . This model was generated using 4 L to 8 L 
data of S. bacillaris strain siva2011 and shows only small 
discrepancies between measured and predicted data. The 
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nondimensional DW* can be converted into dimensional 
DW in g/L by multiplying by the maximum DW experi- 
mentally obtained from the baseline test which in this case 
is 4 L. The nonlinear modeling agrees with measured data 
both qualitatively and quantitatively, where the modeling 
has enhanced the R 2 value up to 95% compared to the lin- 
ear model value 87.4% [16]. This suggests that the nonlin- 
ear regression approach enhances accuracy of modeling 
which provides key scale-up informations of S. bacillaris 
strain siva2011 biomass for RRR-a- tocopherol production. 
In this study, the significance of this empirical approach 
provides insights into the application of a nonlinear re- 
gression model that increases the R 2 value and enhances 
the quantitative predication of S. bacillaris strain siva2011 
scaled-up in a 20 L bioreactor. This allows for the sys- 
tematic understanding and design of a multi-variable 
nonlinear regression experiment for significant biomass 
production. 

Conclusions 

Photosynthetic microalgae are rich in ^ftR-a-tocopherol 
and a potential source for this natural antioxidant which 
is an essential human micronutrient. A significant ad- 
vantage of this natural source is the maintenance of the 
specific bioactive form needed for nutrition and the 
elimination of possible issues with potential unknown 
or unexpected toxicities from the synthetic conforma- 
tions. S. bacillaris strain siva2011 has a unique vitamin 
E biosynthetic mechanism capable of sustaining high 
levels of production, including inducible enhanced pro- 
duction which could provide a possible production plat- 
form of i?^-a-tocopherol for pharmaceutical industries. 
A nonlinear mathematical model was developed to model 
scale-up to production in 20 L reactors using an approach 
with a higher accuracy based on the dual variables tested 
in 4 L and 8 L reactors. The R 2 value from this study dem- 
onstrates that this nonlinear approach significantly im- 
proves estimation of S. bacillaris strain siva2011 biomass 
production in the bioreactor than does the linear model. 
Additional studies with progressively larger reactors (and 
models) will be needed to bridge the gap between labora- 
tory and industrial scale. Nevertheless, this data provides 
enhanced bioprocess engineering information in the pro- 
gression towards large-scale pharmaceutical .ftftR-a-tocoph- 
erol production from S. bacillaris strain siva2011 biomass. 

Methods 

Bioreactor culture, elicitation and analytics 

S. bacillaris strain siva2011 cells were cultured in a balloon 
type bioreactor (4 L and 8 L). The bioreactor experimen- 
tal design and the biomass harvesting were performed 
as described by Sivakumar et al. [16]. For elicitation, 
three concentrations of filter sterilized Meja 50, 100, 
and 200 (iL/L were added to the S. bacillaris strain 



siva2011 culture on the 3 rd day. Meja dissolved in etha- 
nol and Meja not dissolved in ethanol were tested, and 
both had a similar effect. Elicitations were carried out 
for 24, 48, or 72 hrs. Algal cells were harvested and 
freeze-dried according to the Sivakumar et al. [16] 
method. One gram of Meja elicited and unelicited 
freeze-dried algal cells were used for analysis of RRR-a- 
tocopherol. ^ftR-a-tocopherol was processed according 
to the Sivakumar et al. [4] method. The reversed phase- 
high performance liquid chromatography chromato- 
gram was acquired using a Dionex 3000 (LPG3400A) 
system (Thermo Scientific, Sunnyvale, CA) equipped 
with a thermostatted UltiMate 3000 autosampler and a 
Dionex RF 2000 fluorescence detector. The system was 
monitored by the softwares chromeleon (6.80) for in- 
strument control and data acquisition, data reproces- 
sing, and solute quantification, respectively. An Agilent 
Zorbax Eclipse XDB-C8 (250 mm x 4.6 mm, mean par- 
ticle size 5 \im) column or CI 8 (250 mm x 4.6 mm, 
mean particle size 5 um) was used to separate RRR-a- 
tocopherol. The mobile phase consists of a linear gradi- 
ent of 90% methanol in water. The flow rate was 1 ml/ 
min. The total acquisition time was 35 min. The wave- 
length was set at 290 nm for excitation and 330 nm for 
emission. The authenticated &£R-a-tocopherol fluores- 
cence spectra and retention time was used for HPLC 
confirmation of i?i?i?-a-tocopherol in the samples. The 
i^i^-a-tocopherol liquid chromatography mass spec- 
trometry spectrum confirmation was acquired using a 
Shimadzu 8050 mass spectrometer. 

Nonlinear regression 

For a stepwise approach, the first phase of the nonlinear 
regression method begins with the simplification from 
five variables to two variables, time (t) and reactor vol- 
ume (V), as shown in equation 1. The regression will be 
extended to multivariables when the two variables 
method is validated. 

DW =f (y^EC, OrP.pH, t, V)~f (f, V) (1) 

The associated variables DW, t, and V in equation 1 
are nondimensionalized into DW*, t*, and V* as shown 
in equation (2) to (4). 



DW 
DWZ 



V[L] 
4L 



(2) 
(3) 
(4) 



Where DW max is the maximum DW [g/L] produced 
from 4 L test under 0.02% C0 2 fraction, and t max is the 
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maximum time to reach the maximum DW. In this 
study, DW max and t max were 3.45 g/L and 6 days, respect- 
ively. The volume was standardized by the baseline: 4 L in 
this modeling. 

A nonlinear model was assumed by combining 4 th 
order polynomials and the power form of the equation as 
shown in equation (5) in consideration of measured data 
of 4 Land 8 L. 

DW* = (at* 4 + bt* 3 + ct* 2 + dt* 2 + et +/) V* g (5) 

All associated constants a, b, c, and d in equation (5) 
were determined from a nonlinear regression method as 
shown in equation (6). 

DW* = (2.686£* 4 -6.942£* 3 + 5.109£* 2 

+ 0.127* + 0.0184) V* 0 - 151 (6) 

Equation (6) allows the predicting of DW in reactor 
volume 20 L at 0.02% C0 2 fraction from the correlation 
obtained from 4 L and 8 L measured data. 

Statistical analysis 

Bioreactor culture, elicitations, and analytical experiments 
were repeated at least three times, each with three replica- 
tions. Results were expressed as the mean with standard 
errors. Stepwise dual nonlinear regressions were used to 
investigate the relationship between 4 L and 8 L in order 
to predict 20 L bioreactor scale-up. 
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